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Introduction:  Amphiphilic  molecules,  including  phospholipid  monolayers,  have  been 
employed  as  models  for  biological  membranes  in  the  study  of  membrane  structural 
transitions  and  membrane  transport  [1-6].  Due  to  the  mechanical  stability  of  these  layers 
and  their  strong  resistance  to  oxidation  or  reduction,  electrochemical  techniques  have 
proven  quite  valuable  in  the  investigation  of  the  electrical  propenies  of  these  model 
membranes.  Of  particular  note  is  the  work  of  Miller  and  co-workers  and  Lecompte  and 
coworkers  [7-11].  These  investigators  studied,  employing  impedance  measurements  and 
cyclic  voltammetry,  the  electrochemical  behavior  of  a  variety  of  molecules  such  as  pro¬ 
thrombin  at  mercury  electrodes  coated  with  amphiphilic  monolayers  and  in  some  cases 
directly  onto  the  mercury  surface  itself.  Their  studies  suggested  that  there  was  an  electrode 
potential- induced  reorientation  of  the  amphiphilic  molecules  adsorbed  onto  mercury. 

In  more  recent  investigations,  voltammetric  techniques  have  been  used  to  study 
potential-induced  structural  changes  of  phospholipid  monolayers  that  are  irreversibly 
adsorbed  onto  a  Hg  electrode  [12-14].  For  example.  Figure  1  presents  a  steady-state 
cyclic  voltammogram  in  a  0.20  M  KCl  solution  (at  a  sweep  rate  of  100  V/sec)  for  a  Hg 
electrode  coated  with  an  irreversibly-adsorbed  monolayer  of  dioleoylphosphatidylcholine 
(DOPC,  see  structure  below). 
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Although  no  redox  active  center  (within  the  potential  range  of  0.0  to  -1.5  V  vs  SSCE)  is 
present  in  these  molecular  assemblies,  two  very  sharp,  and  reversible  waves  are  observed 
with  peak  potentials  of  -0.90  and  -1.1  V,  respectively. 

The  peak-shape  voltammetric  response  shown  in  Fig.  1  is  typical  of  Hg  electrodes 
coated  with  electroinactive  long-chain  carboxylic  acids  [15],  as  well  as  other  phospholipids 
[12,13].  The  number  of  voltammetric  waves  observed  and  the  electrode  potential  at  which 
they  occur  are  specific  to  the  chemical  nature  of  the  monolayer.  However,  in  all  of  these 
cases,  the  voltammetric  response  is  unquestionably  associated  with  a  conformational 
change  within  the  monolayer,  since  the  molecules  that  comprise  the  monolayer  are 
electroinactive  within  the  investigated  potential  range.  We  find  evidence  for  such  behavior 
in  a  variety  of  publications  where  amphiphilic  redox-active  molecules  adsorbed  onto 
mercury  surfaces  exhibit  additional  voltammetric  features  that  are  not  associated  with  the 
redox  center  in  question.  Recent  examples  of  these  include  the  work  of  Schiffrin  [16]  and 
co-workers  on  the  electrochemical  behavior  ubiquinone  (UQio)  adsorbed  onto  mercury 
electrodes  as  well  as  the  work  of  Tokuda  and  co-workers  [17,18]  and  Camacho  and  co¬ 
workers  [19]  on  long  chain  viologens.  In  all  of  these  cases  there  is  a  voltammetric  feature, 
not  associated  with  the  redox  active  center,  which  can  be  ascribed  to  the  afore-mentioned 
potential-induced  conformational  transitions. 

Although  there  is  a  qualitative  understanding  that  a  transient  current  will  result 
whenever  there  is  a  potential-induced  re-distribution  of  surface  charge  (the  above  response 
for  adsorbed  DOPC  (Figure  1)  clearly  illustrates  such  behavior),  to  our  knowledge,  no 
theoretical  treatment  exists  which  describes  the  current-voltage  response  associated  with 
such  a  process.  In  the  specific  case  of  phospholipid  monolayers  on  Hg  electrodes,  there 
have  been  limited  attempts  at  developing  a  theoretical  framework  that  allows  interpretation 
of  the  voltammetric  response.  Nelson  and  Leermakers  used  a  statistical  model  based  on 
Flory-Huggins  (F-H)  interaction  parameters  to  elucidate  the  potential-dependent  structures 
of  phospholipid  monolayers  and  bilayers  on  Hg  in  equilibrium  with  a  solution  containing 


3 


the  phospholipid  molecule  [13].  Based  on  their  computations,  these  authors  suggested  that 
the  first  voltammetric  peak  in  Figure  1  “corresponds  to  a  competition  between  head  groups 
and  hydrocarbon  tails  for  access  to  the  interface  giving  rise  to  an  inhomogeneous  layer  of 
two  phases  consisting  of  a  thin  bilayer  and  thin  monolayer  respectively”  [13b],  consistent 
with  the  reorientation  of  some  fraction  of  the  molecules  comprising  the  film  from  a  tails- 
down  to  head-group-down  conformation.  This  behavior  is  rationalized  by  assuming  that 
the  affinities  of  the  solvent  (H2O),  head  group,  and  apolar  hydrocarbon  tails  for  the  surface 
are  strong  functions  of  the  surface  hydrophobicity,  which  is  electrode  potential  dependent 
and  decreases  as  the  electrode  potential  is  moved  away  from  the  potential  of  zero  charge 
(Epzc).  Although  this  description  is  quite  reasonable,  unfortunately,  no  structurally- 
sensitive  experimental  measurements  (e.g.,  in-situ  STM  or  x-ray  diffraction  methods)  have 
been  made  that  either  support  or  refute  the  proposed  potential-dependent  structures.  In 
addition,  values  of  the  F-H  interaction  parameters  for  the  various  potential-dependent 
surface  interactions  are  not  known  [13a]. 

More  recent  theoretical  and  experimental  investigations  by  Wangnerud  and  Jonsson 
have  focused  of  the  physical  nature  of  interactions  that  give  rise  to  adsorption  of 
amphiphiles  as  bilayers  on  charged  surfaces,  including,  among  several  factors,  the 
influence  of  electrostatic  interactions  between  the  surface  and  charged  adsorbed  molecules 
[20].  In  this  work,  electrostatic  effects  are  explicitly  treated  using  the  Poisson-Boltzmann 
equation,  yielding  testable  predictions  concerning  the  nature  of  adsorbed  structures  as  a 
function  of  the  surface  charge  density,  electrolyte  concentration,  and  valence  of  the 
amphiphile  co-ion.  A  key  conclusion  of  the  Wangnerud  and  Jonsson  investigations, 
supported  by  in-situ  ellipsometry,  is  that  the  structure  of  the  amphiphile  layer  is  strongly 
dependent  on  electrostatic  effects,  tending  to  form  bilayers  at  high  surface  charge  densities. 

Following  the  general  approach  of  Wangnerud  and  Jonsson,  a  peak-shaped  current 
in  the  voltammetric  response  of  amphiphilic  monolayers  irreversibly-adsorbed  onto  Hg  can 
be  anticipated  by  considering  the  electrostatic  forces  acting  between  the  charged  metal 
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surface  and  the  ionized  head  groups  on  a  lipid  layer.  Consider,  for  instance,  the  idealized 
structures  shown  in  Figure  2  of  a  lipid  monolayer  in  contact  with  a  metal  electrode  and 
immersed  in  an  electrolytic  solution.  For  simplicity,  we  assume  that  the  molecule  has  a 
positively-charged  head  group  and  that  it  is  possible  for  the  molecule  to  rapidly  change 
between  orientations  in  which  the  head  group  is  at  the  metal/lipid  interface  or  at  the 
lipid/solution  interface.  When  the  potential  of  the  electrode  is  sufficiently  positive  of  the 
potential  of  zero  charge  (E  >  Epzc).  such  that  the  electrode  surface  has  a  net  positive 
charge  density,  there  will  be  a  tendency  for  the  molecules  to  orient  such  that  the  head 
groups  are  located  at  the  lipid/solution  interface,  thereby  reducing  the  electrostatic  repulsion 
between  the  surface  and  lipid  layer.  Conversely,  when  the  electrode  potential  is  negative  of 
the  Epzc.  the  electrostatic  attraction  between  the  negatively  charged  electrode  and  positively 
charged  molecule  will  tend  to  orient  the  molecule  with  the  head  group  at  the  metal/lipid 
layer  interface. 

The  reorientation  of  the  lipid  layer  results  in  the  movement  of  a  net  charge  through 
the  electric  field  that  exists  across  the  lipid  layer.  For  the  phospholipids  corresponding  to 
the  data  in  Figure  1,  the  effective  distance  that  charge  is  moved  is  on  the  order  of  the  length 
of  the  molecule  (~25A).  This  displacement  of  charge  will  alter  the  preexisting  electric  field, 
resulting  in  a  measurable  flow  of  current  at  the  electrode  surface.  This  current  is  associated 
only  with  charging  of  the  interface  required  to  minimize  the  electrical  energy  (and,  thus, 
overall  free  energy)  and  does  not  involve  any  redox  chemistry. 

Although  the  structures  and  conformational  transition  depicted  in  Figure  2  are 
highly  speculative,  and  perhaps  unrealistic  based  on  energetic  considerations,  we  will 
employ  this  model  for  the  sole  purpose  of  demonstrating  a  general  method  of  computing 
the  voltammetric  current  resulting  from  a  potential-dependent  conformational  change.  As 
will  be  discussed  in  a  later  section,  the  methodology  described  here  may  be  readily  adapted 
to  other  structures  (including  amphiphilic  bilayers),  as  well  as  other  potential-dependent 
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interactions  (such  as  solvent  dipole  reorientation),  allowing  the  voltammetric  response  to  be 
computed  for  structures  that  are  deduced  from  future  experimental  measurements. 

As  noted  above,  no  theoretical  treatment  exists  which  allows  prediction  or 
interpretation  of  the  voltammetric  response  associated  with  such  a  conformational 
transition.  Consequently,  prior  voltammetric  investigations  of  amphiphilic  monolayers 
have  not  provided  quantitative  information  concerning  the  structural  features  of  the 
monolayer  nor  the  dynamics  of  potential-dependent  conformational  changes.  On  the  other 
hand,  voltammetric  techniques  for  investigating  electroactive  monolayers,  e.g.,  self- 
assembled  monolayers  with  pendant  redox  groups  [21],  are  well-established  and  provide  a 
quantitative  means  of  measuring  kinetic  and  thermodynamic  parameters  associated  with 
electron-transfer,  the  number  density  of  electroactive  molecules  within  the  film,  and  the 
chemical  stability  of  the  monolayer  in  different  oxidation  states  [22].  In  addition,  recent 
theoretical  and  experimental  investigations  have  shown  that  the  voltammetric  response  of  an 
irreversibly-adsorbed  electroactive  monolayer  [23],  and  monolayers  containing  acid/base 
groups  [24],  are  very  sensitive  to  structural  details  of  the  adsorbed  molecule  (e.g.,  chain 
length,  valence)  as  well  as  to  parameters  that  influence  the  interfacial  potential  distribution 
(e.g.,  ionic  strength  of  the  contacting  solution).  These  latter  developments  suggest  that 
similar  structural  information  concerning  amphiphilic  films  might  be  obtained  by 
electrochemical  measurements. 

Based  on  the  success  of  voltammetric  techniques  in  characterizing  electroactive 
films,  our  laboratories  have  recently  initiated  an  investigation  of  the  factors  that  determine 
the  pieak-shaped  voltammetric  response  of  nonelectroactive  monolayers,  with  the  prospect 
of  applying  the  results  to  amphiphilic  monolayers.  This  report  describes  a  first  approach 
based  on  integrating  electrostatic  and  thermodynamic  arguments  describing  the  film 
conformations,  with  the  relevant  equations  describing  the  perturbation  and  response 
functions  of  the  voltammetric  experiment.  Because  of  the  mathematical  complexities  that 
underlie  an  accurate  description  of  the  film  structure,  especially  at  the  microscopic  level,  as 
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well  as  those  associated  with  the  electroanalytical  method,  a  number  of  simplifications  have 
been  made  in  order  to  yield  a  tractable  problem.  Nevertheless,  the  results  represent  the  first 
quantitative  model  of  the  voltammetric  response  of  an  electroinactive  amphiphilic 
monolayer,  and  are  conceptually  useful  in  understanding  how  the  peak- shaped  current 
response  is  related  to  the  structure  of  the  monolayer  and  its  electrostatic  environment. 

Experimental.  Dioleoylphosphatidylcholine  (DOPC,  approximately  99%)  was  used  as 
received  from  Fluka.  Solutions  of  0.20  mg  of  DOPC  per  mL  of  chloroform  were  prepared 
just  prior  to  use.  Electrolyte  solutions  were  prepared  with  ultrapure  (at  least  99.99+%) 
salts  (Aldrich),  which  were  dissolved  (0.20  M)  in  water  which  was  purified  with  a  Milli-Q 
system  and  buffered  with  high  purity  phosphate  salts  or  hydrochloric  acid  (all  from 
Aldrich)  depending  on  the  desired  pH. 

Electrochemical  experiments  were  carried  out  with  an  EG&G  PARC  173 
potentiastat  and  175  Universal  Programmer.  Data  were  collected  on  a  Nicolet  4094  digital 
oscilloscope  and  transferred  to  a  personal  computer  for  analysis. 

A  sodium-saturated  calomel  electrode  (SSCE)  was  used  as  the  reference  electrode 
and  a  large  area  platinum  coil  was  used  as  the  counter  electrode.  The  working  electrode 
was  a  Kemula  type  (micrometer  driven)  hanging  mercury  drop  electrode  (HMDE),  which 
was  filled  with  high  purity  mercury  (electronic  grade,  99.9998%  from  Johnson  Matthey). 
A  two-compartment  cell  without  frits  was  employed  in  all  experiments.  The  cell  had 
standard  joints  for  all  electrodes  and  for  degassing. 

The  detailed  experimental  procedures  have  been  described  previously  [12,  14]. 
Figure  1  shows  the  typical  voltammetric  response  (at  100  V/sec)  for  a  DOPC  monolayer 
adsorbed  on  a  Hg  electrode  surface  in  a  0.20  M  KCl  solution  (pH=8.2).  Two  very  sharp 
and  reversible  waves  centered  at  about  -0.9  and  -1.0  V,  respectively,  are  clearly  evident.  A 
qualitatively  similar  response  was  observed  over  the  pH  range  from  0.7  to  9.2.  However, 
the  potentials  at  which  these  two  waves  appeared  as  well  as  their  sharpness  were  dependent 
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on  pH  and  on  the  nature  of  the  electrolyte.  Such  effects  and  trends  have  been  reported 
earlier  [12]. 

Results  and  Discussion.  An  expression  that  describes  the  voltammetric  behavior  of 
electroinactive  lipid  layers  is  developed  in  the  following  order.  In  Section  I,  an  electrostatic 
model  of  the  interface  that  incorporates  the  essential  features  of  the  assembled  phospholipid 
layer  is  introduced.  General  analytic  expressions  are  presented  which  relate  the  applied 
electrode  potential  in  a  voltammetric  experiment  to  the  orientation  of  the  lipid  molecules  on 
the  surface  and  to  other  electrostatic  parameters  that  influence  the  overall  interfacial  charge 
and  potential  distributions  (e.g.,  electrolyte  concentration  and  dielectric  constant).  In 
Section  II,  an  expression  for  the  free  energy  of  the  lipid  molecules  is  introduced  which 
allows  the  configuration  of  the  molecular  layer  to  be  computed  as  a  function  of  the 
interfacial  potential  distribution.  The  electrostatic  and  thermodynamic  expressions  are 
combined  in  Section  III  to  obtain  the  general  current-potential  response  expected  in  a 
voltammetric  experiment.  In  Section  IV,  theoretical  voltammograms  are  presented  for  a 
variety  of  system  parameters  (e.g.,  surface  coverage,  molecular  dimensions,  degree  of  lipid 
ionization),  and,  in  Section  V,  the  results  are  compared  to  experimental  data. 

The  general  strategy  used  in  developing  a  model  for  the  voltammetric  response  of 
these  lipid  monolayers  is  similar  to  that  reported  earlier  for  analyzing  the  behavior  of  the 
monolayers  containing  redox  or  acid/base  functionalities  [23,24], 

I,  Inteifacial  Potential  and  Charge  Distribution.  Figure  3  shows  the  electrostatic  model 
used  to  compute  the  interfacial  potential  distribution.  The  model  consists  of  a  tightly- 
packed  lipid  monolayer  assembled  on  a  planar  metal  electrode  of  area  A  (cm^)  and 
immersed  in  an  electrolyte  solution  that  contains  a  symmetric  z:z  electrolyte.  The  solution 
is  assigned  a  dielectric  constant, 
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The  lipid  monolayer  assembly  can  be  modeled  in  several  ways  that  strongly  depend 
on  the  assumed  molecular  conformation  and  the  degree  of  dissociation  of  head  group 
functionalities.  For  instance,  at  moderate  and  high  pH,  DOPC  exists  as  a  zwitterionic 
species  (see  structure);  whereas  at  low  pH  (pH<2),  the  phosphate  group  is  protonated, 
leaving  the  head  group  with  a  net  positive  charge.  As  mentioned  earlier,  the  voltammetric 
response  is  slightly  dependent  on  pH,  but  two  reversible  voltammetric  peaks,  similar  to  the 
response  shown  in  Figure  1,  are  observed  at  all  pHs. 

The  model  we  assume,  and  schematically  shown  in  Figure  3,  is  that  of  an  electrode 
surface  covered  by  two  dielectric  films  of  thicknesses  di  and  dz  with  dielectric  constants, 
El  and  £2,  respectively.  The  distance  di  is  defined  as  the  distance  between  the  surface  and 
the  positive  site  when  the  lipid  head  group  is  located  at  the  electrode/lipid  interface.  It  is 
difficult  to  obtain  a  precise  value  for  di  without  detailed  structural  information  concerning 
the  orientation  of  the  head  group  at  the  surface.  However,  based  on  the  molecular  structure 
shown  earlier,  it  is  likely  that  this  distance  has  a  value  between  2  and  3  A.  The  thickness  of 
the  dielectric  layer  d2  is  defined  as  the  length  of  the  molecule,  minus  the  length  di.  A 
transition  between  the  inwards  (head  groups  down)  and  outwards  (tails  down)  orientations 
will  result  in  a  displacement  of  the  positive  charge  on  the  molecule  by  a  distance  of  d2  -  di. 

The  surface  coverage  of  lipid  molecules  oriented  with  their  polar  head  groups  at  the 
electrode/lipid  and  lipid/solution  interfaces  are  designated  as  Fi  and  r2,  respectively.  Since 

the  molecules  are  irreversibly  adsorbed,  the  total  surface  coverage.  Ft,  is  defined  asFT  = 

n+Fz. 

To  simplify  the  resulting  expressions,  we  will  assume  each  lipid  molecule 
possesses  a  single  net  positive  charge,  Zmol.  For  the  low  pH  DOPC  structure,  we 
anticipate  Zmol  ~  +1.  From  electrostatic  considerations,  at  high  pH,  the  zwitterionic  form 
of  DOPC  would  be  expected  to  behave  equivalently  to  a  positively  charged  species,  with 
Zmol  loss  than  +1.  This  conclusion  is  based  on  the  dipole  nature  of  the  zwitterion,  in 
which  the  positively  charged  choline  group  undergoes  a  slightly  larger  spatial  translation 
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than  the  negatively-charged  phosphate  group  (~4  -  6  A  based  on  the  molecular  dimensions 
discussed  above)  during  the  assumed  conformational  transition.  Thus,  the  image  charge 
(voltanunetric  current),  induced  in  the  metal  during  the  conformational  transition  of  the 
zwitterion  will  be  equivalent  to  that  generated  by  the  transition  of  a  molecule  containing  a 
small  positive  charge.  However,  these  general  arguments  could  be  somewhat 
compromised  by  the  fact  that  ion  association  between  the  amphiphilic  molecules  and 
dissolved  counter-ions  may  be  significant,  albeit  unknown  [25].  As  detailed  below,  this 
latter  limitation  prevents  an  independent  measurement  of  Zmoi  and  rT  using  voltammetric 
data.  However,  the  charge  surface  density,  Zmoirx,  involved  in  the  conformational 
transition  is  directly  measurable. 

Assuming  that  there  is  no  charge  in  the  film  except  that  associated  with  the  charged 
head  groups  (0  <  Zmol  <  +1),  the  potential  will  decay  linearly  from  (|)ni  to  (l)i  and  from  <l)i  to 
<{>2.  In  the  solution  phase  (x  >  di  +  d2),  the  potential  is  assumed  to  be  governed  by  the 
Gouy-Chapman  model  and  decays  to  (t)s  in  the  bulk  solution: 

tanh(ze(0U)-0,)/2kr)  =  tanh(2e(02 -0,)/2/:r)exp(-K:U-(d, -i-dj))) 

where  <J)(x)  is  the  potential  (V)  in  the  solution  at  position  x,  k  is  the  inverse  Debye  length 
(m-1)  given  by  K  =  [z^e'^2n°/e,£,kTf ,  e  is  the  charge  of  an  electron  (C),  nO  is  the 

number  concentration  of  the  ions  in  the  electrolyte  (cm’^),  Eq  is  the  permitivity  of  free  space 
(C^  cm‘2)  and  k  is  Boltzmann's  constant  (JK‘l). 

In  addition  to  the  above  assumptions,  we  also  assume  that  charges  on  the  metal 
surface  (at  x  =  0)  and  on  the  head  groups  (at  di  and  di+d2)  are  delocalized  (i.e., 
discreteness  of  charge  is  ignored).  Inclusion  of  the  finite  size  of  ions  in  the  description  of 
the  electrochemical  behavior  of  redox  and  acid/base  functionalized  monolayers  on 
electrodes  has  been  recently  discussed  by  Fawcett  and  coworkers  [23b, 24b].  As  will  be 
shown  below,  the  potential  distribution  of  the  system  we  are  considering  is  dominated  by 
the  capacitance  of  the  dielectric  region  corresponding  to  the  hydrocarbon  tail  (i.e.,  the 
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region  defined  by  d2  -  di).  Thus,  the  approximation  of  delocalized  charges  outside  of  this 

region  (i.e.,  at  the  metal  and  solution  interfaces)  appears  reasonable. 

The  electrostatic  force  exerted  on  the  lipid  molecules  at  any  applied  electrode 

potential  can  be  obtained  by  considering  the  relationship  between  the  potential  distribution 

and  the  charge  densities  on  the  metal,  at  positions  di  and  (di  +  d2),  and  in  the  electrolyte 

solution.  Since  the  electrode  surface  area  is  generally  of  the  order  of  cm^,  and  the 

interfacial  distances  di,  d2,  and  k-i  are  of  the  order  of  nm,  it  is  only  necessary  to  consider 

the  potential  distribution  normal  to  the  surface.  With  attention  drawn  to  Figure  3  the 
dV 

electric  field  (E  =  -  within  each  region  of  the  interface  is  given  by  eq.l; 


E  =  0  (*.•  =  oo) 

ze  L  IkT 

0 

The  electric  fields  and  the  charge  densities  associated  with  the  metal,  lipid  layer,  and 
solution  are  related  through  Gauss'  law.  The  Gaussian  boxes  shown  in  Fig.  3  enclose  the 
charges  of  interest.  For  the  assumed  planar  geometry,  the  electric  field  passes  only  through 
the  faces  of  the  boxes,  each  of  surface  area  S.  Thus: 


x  <0 
0<x<d, 
di  <  jc  +  ^4) 
x>d^  +^2 

ac^oo  (1) 


Combining  eq.  1  and  2  and  applying  the  result  to  the  four  boxes  shown  in  Figure  3,  yields 
the  following  system  of  equations  relating  the  charge  densities  (Oi  =  Qi/S)  and  interfacial 
potentials. 


box  1 

(3a) 

box  2 

(3b) 

box  3 

(3c) 

O'm  +  O-j  +  0-2  -1-  Ojif  =  0 

box  4 

(3d) 

where  the  definitions  of  charge  densities  (C‘cm*^)  on  the  electrode  surface  (<Jrn)>  at  the 
planes  of  the  lipid  head  groups  (ai,  02),  and  in  the  solution  (the  diffuse  layer)  (adif)  are 
given  by  eq.  4 


x  =  0 

X  =  d 

02 

X  =  di  +  d2 

Odif 

di  -t-  d2  <  X  <  00. 

As  has  been  indicated  earlier,  the  lipid  layer  is  electrochemically  inactive  and  irreversibly 
adsorbed  over  the  potential  range  of  interest  Thus,  the  total  charge  per  unit  area  in  the  film 
is  constant,  i.e.. 


+  =  (5) 

n.  Relationship  Between  Local  Potential  and  Molecular  Orientation.  In  order  to  calculate 
the  voltammetric  response  corresponding  to  the  field-induced  transition,  it  is  necessary  to 
establish  a  relationship  between  the  two  assumed  molecular  orientations  and  the  electrode 
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potential.  As  previously  noted,  the  molecular  orientation  is  anticipated  to  be  a  function  of 
the  applied  potential  if  the  electrostatic  forces  between  the  charges  on  the  metal  surface  and 
on  the  lipid  molecule  are  sufficiently  large  to  flip  the  molecules  between  the  inwards  and 
outwards-directed  configurations,  Figure  2.  The  fact  that  a  reversible  voltammetric 
response.  Figure  1,  is  observed  even  at  relatively  high  scan  rates  (e.g.,  100  V/s)  suggests 
that  this  process  is  facile  on  the  timescale  of  the  electrochemical  experiment,  and  that  the 
equilibrium  structure  is  determined,  to  a  significant  extent,  by  the  interfacial  potential 
distribution. 

To  account  for  the  influence  of  the  electrostatic  potential  on  the  molecular 
orientation,  we  assume  that  the  electrochemical  potential  of  the  lipid  molecules  is  given  by 

Jij  =  // J  -f-  In  aj  +  z^iF<l>j  (6) 

where  F,  R  and  T  are  Faraday's  constant  (C  moF^),  the  molar  gas  constant  (J  moF^  K'l), 
and  the  absolute  temperature  (K),  respectively.  The  subscript  j  refers  to  the  two  possible 
molecular  orientations  (inwards  (j  =  1)  and  outwards  (j  =  2))  of  the  molecules  adsorbed  on 
the  electrode  surface.  The  constants  /tf  and  /tf  indicate  non-electrostatic  energy  terms 

corresponding  to  hydrophobic/hydrophilic  interactions,  dispersion  and  solvation  forces, 
and  steric  effects  associated  with  the  amphiphilic  layer  in  the  inwards  and  outwards 
orientations,  respectively.  We  treat  these  terms  as  being  independent  of  the  electrode 
potential,  since  the  predominant  potential-dependent  interactions  with  the  electrode  surface 
are  likely  to  involve  electrostatic  forces  between  the  surface  charge  and  the  fixed  charge  on 
the  head  group.  The  reader  is  referred  to  the  thorough  discussion  of  Wangnerud  and 
Jonsson  [20]  for  a  detailed  physical  description  and  estimation  of  the  non-electrostatic  work 
terms. 
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At  equilibrium,  /i,  =  /i2.  Approximating  the  activities  ai  and  a2  by  their  respective 
surface  concentrations,  Fj  and  r2,  yields  the  relationship  between  the  electrostatic 
potential  distribution  in  the  film  and  the  orientation  of  the  charged  sites: 


01  02  ~ 


^molF 


InlL. 

r. 


^molF 


(7) 


Eq.  (7)  is  simplified  by  defining  /  =  fj  / r^.,  and  the  potential, E"  =  -(/xf  -^1)1  z^^F, 
yielding 


0,-02  =  ^"— ^In 


^  f  ^ 

,1-/ 


(8) 


It  is  clear  that  the  potential  E®  represents  the  difference  in  the  standard- state  chemical 
potentials  of  the  assumed  molecular  conformations,  and,  in  this  sense,  its  use  is  similar  to 
that  employed  in  defining  the  standard  redox  potentials  of  electron-transfer  reactions  (where 
11°  and  /if  represent  the  standard-state  chemical  potentials  of  the  reduced  and  oxidized 
partners  of  the  redox  couple).  As  detailed  in  a  later  section,  the  peak  position  of  the 
voltammetric  wave  corresponding  to  a  conformational  transition  is  primarily  determined  by 
the  value  of  Eo  (versus  a  yet  specified  reference  potential).  Experimental  variables  that 
affect  the  electrostatic  potential  distribution  (i.e.,  dielectric  constant,  ion  concentration)  have 
a  smaller,  but  significant  effect,  on  the  potential  at  which  the  structural  transition  is 
observed. 

Admittedly,  the  model  that  we  have  assumed  for  the  potential-dependent 
conformational  transition  may  not  be  necessarily  correct,  and  clearly  does  not  include  all  of 
the  features  necessary  to  describe  the  actual  amphiphilic  layer.  However,  although  it  is 
beyond  the  scope  of  this  preliminary  treatment,  both  the  electrostatic  and  thermodynamic 
treatments  outlined  above  may  be  readily  modified  when  considering  different  structures 
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and  potential  dependent  interactions.  For  instance,  the  Gaussian  box  approach  used  above 
can  be  readily  extended  to  include  a  true  zwitterionic  structure,  requiring  only  the  addition 
of  an  additional  charged  layer  in  the  model  schematic  (Figure  3)  and  one  additional 
application  of  Gauss’s  law  (eq.  (2)).  In  analogous  fashion,  the  electrostatic  model  can  be 
extended  to  bilayer  structures,  and  ion-association,  and  finite  ion  size  effects  may  be 
included,  as  was  recently  done  by  Fawcett  for  self-assembled  redox  and  acid/base 
monolayers.  Similarly,  the  thermodynamic  expressions  we  assume  are  highly  idealized  in 
many  respects.  A  key  deficiency  is  the  implicit  assumption  that  the  energetics  of 
reorientation  of  any  molecule  within  the  film  are  independent  of  the  state  of  orientation  of 
the  neighboring  molecules,  i.e.,  the  probability  of  a  molecule  being  in  an  up  or  down 
configuration  is  independent  of  the  orientational  state  of  the  neighboring  molecules.  Lateral 
electrostatic  interactions  within  the  film  will  likely  make  the  random  approximation 
improbable.  This  problem  can  be  treated  by  methods  similar  to  those  developed  by  Levine 
[26]  for  accounting  for  lateral  dipole-dipole  interactions  within  adsorbed  solvent  layers,  or 
in  a  more  heuristic  fashion  by  employing  an  adjustable  interaction  parameter  in  eq.  (6) 
(equivalent  to  employing  a  Frumkin-type  isotherm). 

in.  Voltammetric  Response. 

If  the  electrode  potential  is  scanned  at  a  constant  rate  v  (V/s),  the  voltammetric 
response  of  an  electrode  coated  with  an  electrochemically-inert  film  is  due  only  to  the 
capacitive  charging  current 

i  =  vACj  (9) 

In  eq.  (9),  Cj  is  the  total  capacitance  (Fcm'^)  of  the  interface,  which  can  be  evaluated 
using  the  formal  definition:  Cj  =  da^  I  dE  =  da„  /  In  deriving  an  expression  for 
Ct,  it  is  convenient  to  define  the  following  interfacial  capacitances: 
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C\  —  ^0^1 ! 

(lOa) 

(10b) 

Cdif  =  £o^s^^osh[2e{<l)^-(ps)/2kT]. 

(10c) 

Cl  and  C2  represent  the  potential-independent  capacitances  of  the  film  and  Cdif  is  the 
potential-dependent  diffuse  layer  capacitance. 

The  total  interfacial  potential  drop  ((t)m  -  ^s)  is  equal  to  the  sum  of  the  individual 
contributions  across  the  interface: 

+  +  (11) 

Substituting  eqs.  3  and  10  into  eq.  1 1  yields: 

<t>m  =  +  Q'ko-^  +  fT,)  -H  02  (12) 

Differentiating  eq.  12  with  respect  to  Om  yields  the  potential-dependent  interfacial 
capacitance  of  the  system: 


+  C 


-1 


(13) 


Expressing  the  quantity  d<p2  /  d(J„  ^  (d<l>2 1  I  (d<j„  I  ,  and  noting  that 
do^/  da^f  =  -1  (from  eqs.  3d  and  5),  and  id(l>2l  d(Jdif)  =  -C^  (from  eqs.  3c  and 
10c)  yields 

Cf‘  =  Cf‘  -I-  Q'^l  +  (14) 
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An  explicit  expression  for  9(7i/3<Jm  in  terms  of  f  is  derived  from  the  definition  of  cjm. 
Substituting  eq.  8  into  eq.  3b  (along  with  the  definition  =  z^,fF  F^)  yields 


^0^2 

4 


£"  + 


RT 


which  is  differentiated  with  respect  to  ai: 


(9(T,  _  dcx!  df 


g.g2  RT  1  .  , 


(15) 


(16) 


Since  0  ^  f  <  1,  the  bracketed  expression  of  eq.l6  is  always  positive,  and  therefore 
<9(Tj  /  do^  must  always  be  negative.  Physically,  this  relationship  indicates  that  an  increase 

in  the  positive  charge  on  the  electrode  must  result  in  a  decrease  in  the  positive  charge  at  x  = 
di,  and  conversely,  that  an  increase  in  the  negative  surface  charge  will  result  in  an  increase 
in  the  positive  charge  at  di. 

Combining  eqs.  9,  14,  and  16  yields  the  voltammetric  current  for  the  amphiphilic 
monolayer-coated  electrode  in  terms  of  f. 


V - 

II 

Cf'  +  C2 

f  p 

1  - 

1 

1 

RT 


^  z 


J 

moi 


/(I-/) 


+  1 


^-1 


(17) 


Before  proceeding  to  the  numerical  calculations,  it  is  instructive  to  note  that  the 
voltammetric  current  will  be  dominated  by  the  smallest  individual  capacitance  of  the 
interface.  Consider,  for  example,  a  lipid  layer  with  di  *  4  A,  d2  =  30  A,  and  ei  =  £2  =  4, 
immersed  in  an  aqueous  solution  containing  0.10  M  of  a  1:1  supporting  electrolyte  (Eg  =  78 
and  K-1  ~10  A).  From  eq.  10,  Q  =  10//F  /  cm^  C2  =  1.5//F  /  crn^  and  -  90/iF  /  cm^ 
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(where  is  the  minimum  diffuse  layer  capacitance  evaluated  at  <|)2  =  (t)s).  Thus,  the 
voltammetric  current  will  be  dominated  by  the  capacitance  associated  with  the  hydrocarbon 
tail,  C2,  and  the  variation  of  +  /  d<y„),  eq.  14,  with  the  electrode  potential  will 

largely  determine  the  shape  of  the  voltammetric  response.  A  maximum  in  the  voltammetric 
curve  is  expected,  since  as  the  film  structure  changes  from  a  fully  inwards  orientation  (f  = 
1)  to  a  fully  outwards  orientation  (f  =  0),  the  quantity  (da^  /  daj  changes  from  0  to  a 
negative  value  (see  eq.  16),  and  +  d<7^  /  d(T„)  decreases. 

As  the  film  becomes  either  completely  head  group  directed  inward  or  outwards 
relative  to  the  electrode  surface  (i.e.,  as  /  ->  1  or  /  -4  0),  d<j^  /  ->  0  and  eq.  14 

reduces  to  the  expression  for  the  total  capacitance  of  a  film  in  which  no  conformational 
changes  are  present. 

Cf  =  C,  +C2’ +C^  (18) 


rv.  Computational  Results. 

The  voltammetric  response  is  obtained  by  calculating  the  current  and  applied 
potential  as  a  function  of  f.  The  current,  i,  is  obtained  directly  from  eq.  (17).  An  explicit 
expression  for  the  electrode  potential,  E  (=  ((Jm  -  ())s),  is  readily  obtained  by  solving  eq.  (3c) 
for  (t)2  (using  the  identity  sinh-l[u]  =  ln[u  +  (u2  +  1)0-5]),  and  substituting  the  result  into  eq. 
12,  along  with  the  definitions  oi  =  ZmolfFTr,  Oi  +  02  =  ZhioiFTt,  and  Odif  =  -(am  + 
ai  +  02)  (eq.  3d).  The  result  is  given  in  eq.  (19),  where  Om  (eq.  (15))  is  computed  for 
each  specified  value  of  f. 


E  =  C-  +  C-\g„  +  z^ifFTr)  + 


(19) 


(2kT/ze)ln 


ze(G„+z^FTr)  1 

0.5^ 

2e^£^KkT 

(  2e^e,KkT  )  _ 

j 
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Eqs.  (17)  and  (19)  completely  describe  the  voltammetric  response  in  terms  of  the  properties 
of  the  lipid  layer  (surface  coverage,  charge,  molecular  dimensions  and  dielectric),  the 
solution  phase  (electrolyte  concentration,  charge,  and  solution  dielectric),  and  temperature. 
The  following  independent  parameters  are  input  into  each  calculation:  Eo,  Eoei/di  (=  Ci), 
eo£2/d2  (=  C2),  ZtnoirT,  z,  T,  V,  and  A.  All  calculations  presented  are  computed  for  T 
=  300  K  and  z  =  1.  Unless  stated  otherwise,  Eo  is  taken  equal  to  Epzc  (i  e.,  E®  =  0). 
Anodic  and  cathodic  currents  are  mirror  images  of  each  other,  since  the  lipid  layer  is 
assumed  to  be  in  equilibrium  with  the  electrode  potential  at  all  times.  Only  the  cathodic 
current  will  be  presented  in  the  figures.  Electrode  potentials,  E,  are  referenced  to  the 
potential  of  zero  charge,  Epzc,  of  the  bare  electrode  (i.e.,  in  the  absence  of  the  adsorbate;  as 
shown  below,  the  adsorption  of  the  lipid  layer  induces  a  significant  shift  in  the  Epzc). 
Voltammetric  currents  are  normalized  to  the  scan  rate  and  electrode  area  (iAA). 

Figure  4  shows  a  voltammogram  calculated  using  the  following  parameters:  Zmol  = 
0.5,  Ft  =  6  X  10-11  mol/cm2,  ej/di  =  10/2,  £2/^2  =  6/33,  Eg  =  78,  and  k-i  =  9.6  A 
(corresponding  to  a  0.10  M  aqueous  solution  of  a  1:1  electrolyte).  These  values  are  chosen 
to  illustrate  the  general  electrochemical  behavior  of  a  lipid  monolayer-coated  Hg  electrode. 
In  later  sections,  the  effects  of  varying  each  parameter  will  be  examined  in  detail. 

Qualitatively,  the  shape  of  the  theoretical  voltammogram  is  in  good  agreement  with 
the  1st  voltammetric  wave  observed  for  a  Hg/DOPC  electrode  in  a  pH  8.2,  0.20  M  KCl 
solution  (Figure  1).  (As  previously  noted,  the  simplified  model  we  have  assumed  does  not 
account  for  the  formation  of  the  final  bilayer  structure,  proposed  as  the  origin  of  the  2nd 
cathodic  wave  [13];  thus,  hereafter,  all  comparisons  between  experiment  and  theory  refer 
only  to  the  1st  wave.)  A  peak-shaped  voltammetric  response  is  observed,  centered  slightly 
negative  of  Epzc-  As  shown  below,  the  potential  corresponding  to  the  current  maximum 
and  the  waveshape  are  dependent  on  Eo,  Zmoi,  Ft,  Eg,  k"!,  and  the  capacitances  of  the 
lipid  monolayer  (£i/di  and  £2/d2). 
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The  shape  of  the  voltammetric  wave,  Figure  4,  is  similar  to  that  expected  for  an 
electrode  coated  with  an  irreversibly  adsorbed  monolayer  of  an  electroactive  species.  It  is 
thus  especially  interesting  that  the  theoretical  voltammogram  is  entirely  due  to  capacitive 
charging  of  the  interface,  and  does  not  involve  the  transfer  of  electrons.  A  physical 
explanation  of  voltammetric  response  is  illustrated  in  Figure  4,  which  shows  the  fraction 
of  molecules  with  polar  head  groups  oriented  towards  the  interface,  f  =  Fi/Tt  (dashed 
line),  plotted  as  a  function  of  the  electrode  potential.  At  positive  potentials,  i.e.,  E  > 
Epzc.  f  is  essentially  equal  to  0,  indicating  that  all  of  the  lipid  molecules  in  the  film  are 
oriented  with  their  head  groups  oriented  away  from  the  interface.  Under  these  conditions, 
the  capacitance  of  the  electrode  is  essentially  constant,  since  the  total  electrode  capacitance, 
Ct’^  =  Cr^  +  C2'^  +  Cdif  ^  eq.(18),  is  dominated  by  the  smallest  individual  capacitance, 
which  corresponds  to  the  potential-independent  C2.  Consequently,  a  constant  "background 
or  baseline"  charging  current  is  observed  at  values  of  E  sufficiently  positive  of  Epzc-  The 
capacitance  calculated  from  this  constant  current  is  ~2  ^F/cm2,  in  reasonable  agreement 
with  the  experimental  value  obtained  from  Figure  1  (2.5  ^F/cm2).  At  potentials  slightly 
positive  of  Epzc.  the  lipid  molecules  that  comprise  the  film  begin  to  flip  from  an  outwards 
orientation  to  an  inwards  orientation,  due  to  the  decrease  in  the  electrostatic  repulsion 
between  the  positive  head  groups  and  the  positively  charged  surface.  The  quantity  f 
increases  monotonically  as  the  potential  is  scanned  negative  of  Epzc.  s^t  sufficiently 
negative  potentials,  f  approaches  1  as  the  electrostatic  attraction  between  the  head  groups 
and  negative  surface  causes  a  complete  inward-directed  orientation  of  the  molecular  film. 

Further  insight  into  the  origin  of  the  voltammetric  current  is  obtained  by  examining 
the  charge  density  on  the  electrode  surface,  Om,  as  a  function  of  the  applied  potential,  E. 
Fig.  5  shows  Om  vs  E  (computed  using  eqs.  (15)  and  (19))  for  the  same  set  of  parameters 
used  in  computing  the  voltammetric  response  in  Fig.  4.  At  potentials  sufficiently  positive 
or  negative  of  the  Epzc  die  slope  of  Om  vs.  E  is  constant,  indicating,  as  before,  that  the 
electrode  capacitance  (Ct  =  3atn/3E)  is  dominated  by  the  potential-independent  lipid-layer 
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capacitance,  C2.  Near  Ep^c,  Cm  changes  rapidly  in  response  to  the  reorientation  of 
adsorbed  molecules.  The  peak-shaped  voltammogram  of  Figure  4  is  the  result  of  the  flow 
of  electrons  to  the  electrode  surface  (but  not  across  the  interface)  in  response  to  the  field- 
induced  movement  of  the  positive  charge  on  the  lipid  molecule  towards  the  surface. 

In  general,  the  adsorption  of  a  charged  species  from  solution  will  induce  a  charge 
on  the  metal  surface,  the  magnitude  of  which  will  depend  on  all  parameters  that  affect  the 
electrostatic  potential  distribution.  For  the  adsorption  of  a  positively  charged  phospholipid 
monolayer,  the  Epzc  shifts  to  a  more  positive  potential.  The  plot  of  Cm  vs  E,  Figure  5, 
allows  a  graphical  means  of  determining  the  shift  in  Epzc  due  to  adsorption  of  the  lipid 
layer.  The  intersection  of  a  horizontal  dashed  line  corresponding  to  Gm  =  0  with  the  Om  vs 
E  curve  defines  the  Epzc  after  adsorption,  which  we  will  hereafter  refer  to  as  Ep^c-  This 
intersection  occurs  at  0.167  V  vs  Epzc-  Thus,  adsorption  of  the  charged  lipid  molecules 
onto  the  Hg  surface  results  in  a  0.167  V  shift  in  the  potential  of  zero  charge  [27]. 

Figure  6  shows  the  potential  distribution  across  the  lipid-coated  electrode  interface 
for  values  of  f  between  0.00034  and  0.99966,  calculated  using  the  same  set  of  system 
parameters  employed  above.  As  anticipated  from  the  previous  results,  most  of  the  potential 
drop  occurs  across  the  region  defined  by  the  long  hydrocarbon  tail  of  the  phospholipid 
(from  di  to  62).  Regardless  of  the  applied  potential,  the  potential  drop  in  the  solution  phase 
is  small  for  electrolyte  concentrations  >  0.1  M,  a  consequence  of  C2  <  Cdif  for  all  values  of 
f.  Figure  6  also  demonstrates  that  the  Ep>zc  changes  upon  adsorption  of  a  positively  charge 
molecule.  At  E  =  0  (which  corresponds  to  E  =  Epzc  of  the  bare  electrode)  the  electric  field 
at  the  surface  is  seen  to  have  a  finite  negative  value,  as  evidence  by  the  non-zero  positive 
slope  of  the  potential  profile  between  x  =  0  and  x  =  d^.  Thus,  it  is  immediately  apparent 
that  a  finite  negative  charge  exists  on  the  electrode  surface  at  this  potential.  The 
corresponding  shift  in  the  Epzc  to  that  results  from  this  induced  charge  may  be 

calculated  as  described  in  the  previous  paragraph. 
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Figures  7,  8,  and  9  show  the  effects  of  varying  the  following  parameters:  E®,  Fj, 
Celec.  ei/dl,  e2/d2.  Varying  £«  to  be  between  +  0.4  V  has  the  effect  of  shifting  the 
voltammetric  wave  on  the  potential  axis,  Figure  7,  by  an  amount  approximately  equal  to 
E®.  The  linear  dependence  of  the  wave  position  on  E®  is  a  consequence  of  the  fact  that 
virtually  all  of  the  interfacial  potential  drop  occurs  in  the  region  between  di  and  d2,  as 
shown  in  Fig.  6  (i.e.,  ((J)i  -  <t)2) «  E,  neglecting  the  small  potential  drops  between  x  =  0  and 
di,  and  in  the  solution  phase).  The  potential  drop,  (j)!  -  fyi,  determines  the  conformation  of 
the  monolayer,  and  eq.  (8)  shows  that  changing  the  value  of  E®  has  the  effect  of  adding  a 
constant  value  (=E®)  to  the  driving  force  that  is  necessary  to  induce  the  conformational 
transition.  The  key  point  here  is  that  an  experimental  measurement  of  wave  position  (vs 
Epzc)  allows  a  direct  estimation  of  the  free  energy  change  associated  with  the 
conformational  transition.  It  is  noteworthy  that  the  voltammetric  peak  position  is  not 
determined  by  the  sign  of  the  charge  associated  with  the  polar  head  group.  Thus, 
amphiphilic  layers  containing  either  negatively  or  positively  charged  head  groups  may 
exhibit  transitions  at  potentials  negative  of  Epzc-  The  primary  factor  determining  the  peak 
potential,  Ep,  is  the  difference  in  the  non-electrostatic  energy  contributions  (/tf  -  /^l)  of 
the  two  molecular  conformations.  In  addition  to  the  shift  in  Ep,  the  shape  of  the 
voltammetric  wave  has  a  finite,  but  very  weak  dependence  on  the  value  of  E®  (the 
difference  in  the  shapes  of  the  three  waves  in  Figure  8  is  of  the  order  of  the  width  or  the 
pen  used  to  draw  the  curves). 

An  increase  in  Fj  results  in  a  proportional  increase  in  the  voltammetric  peak  height. 
Figure  8a  (measured  relative  to  the  flat  baseline  capacitive  current),  in  agreement  with  the 
predicted  dependence  of  i  on  Ft  indicated  by  eq.  (17).  Increasing  the  charge  on  the  lipid 
molecule,  Zmoi,  has  a  similar  effect  on  the  voltammetric  response  (not  shown)  as  increasing 
Fx-  In  addition  to  increasing  the  peak  height,  an  increase  in  Zfnol  also  tends  to  cause  a 
small  broadening  of  the  wave.  However,  since  the  dependencies  of  the  voltammetric 
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response  on  Zmol  and  Tj  are  so  similar,  there  is  no  simple  means  of  separating  the  effects 
of  these  two  factors. 

The  effect  of  varying  the  supporting  electrolyte  concentration,  Ceiec.  on  the 
voltammetric  wave  is  shown  Figure  8b.  Lowering  Cgiec  induces  a  significant  positive  shift 
in  the  voltammetric  wave,  indicating  that  the  driving  force  for  the  conformational  transition 
is  slightly  larger  at  more  positive  potentials.  This  behavior  can  be  rationalized  by  noting 
that  a  decrease  in  Cgiec  causes  a  larger  fraction  of  the  total  interfacial  potential  drop  to  occur 
across  the  diffuse  double  layer.  From  eq.  11,  an  increase  in  the  diffuse  double  layer 
potential  ({1>2  -  (l)s),  for  any  fixed  electrode  potential  (E  =  (<{)m  -  (^s))  results  in  a  decrease  in 
(01  -  02).  thereby  increasing  the  driving  force  for  the  lipid  molecules  to  reorient  with  their 
head  groups  towards  the  electrode  surface  (see  eq.  (7)  and  Figure  6).  Thus,  there  is  a 
positive  shift  in  the  position  of  the  voltammetric  wave. 

As  shown  in  Figure  9,  the  voltammetric  response  is  a  relatively  strong  function  of 
the  molecular  parameters,  ei/di  and  £2/d2.  For  example,  decreasing  ei/di  results  in  a 
dramatic  decrease  in  the  peak  height,  a  broadening  of  the  wave,  and  a  shift  in  Ep  to  negative 
potentials.  The  dependence  is  analogous  to  that  observed  upon  lowering  Ceiec.  and  is  a 
consequence  of  the  capacitance,  Ci,  (and,  therefore  the  potential  drop,  -  0i)  in  the 
region  between  x  =  0  and  di  becoming  significantly  larger  as  £2/d2  decreases.  Figure  9b 
shows  that  the  primary  effect  of  decreasing  £2/d2  is  to  increase  the  baseline  capacitive 
current  (that  is,  the  effects  in  the  figure  are  real)  defined  by  eq.  (18).  For  sufficiently  large 
values  of  £2/d2,  the  voltammetric  wave  becomes  noticeably  broadened,  a  result  of  C2 
becoming  comparable  to  Ci  and  Cdif. 

For  redox-active  monolayers,  the  coverage  (rx(mol/cm2)  of  electroactive 
molecules  can  be  readily  obtained  from  integration  of  the  charge  defined  by  the 
voltammetric  wave  [22].  No  similar  method  of  analysis  exists  in  the  literature  for 
electrochemically-inert  monolayers.  To  investigate  the  possibility  of  employing  the 
voltammetric  response  in  determining  the  surface  coverage  of  DOPC  on  Hg,  we  have 
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numerically  evaluated  the  charge  under  the  theoretical  voltammetric  waves,  Qapp,  as  a 
function  of  the  charge  corresponding  to  the  actual  coverage  employed  in  the  calculation,  Qf 
(=  ZhioiFTt)  [28].  The  results  are  shown  in  Table  I,  as  a  function  of  the  assumed  values  of 
the  molecular  parameters,  £i/di  and  £2/d2  (other  parameters  used  in  the  calculation  are  the 
same  as  in  Fig.4).  As  seen  in  Table  I,  the  charge  measured  underneath  the  voltammetric 
wave  is  always  smaller  than  the  actual  value  corresponding  to  the  true  surface  coverage. 
However,  the  numerical  examples  presented  in  Table  I  show  that  Qapp  is  never  more  than 
-25%  smaller  than  Qf  over  a  wide  range  of  £i/di  and  £2/d2  combinations.  Thus,  it  is 
possible  to  obtain  reasonably  good  estimates  of  Fj  from  the  voltammetric  response. 

V.  Comparison  of  Experimental  and  Theoretical  Voltammetric  Response  ofHgIDOPC. 

In  Figure  10,  data  from  Figure  1  are  quantitatively  compared  to  a  theoretical 
voltammogram.  Parameters  used  in  computing  the  theoretical  voltammogram  were 
obtained  as  follows.  A  value  of  £2/d2=  0.241  A'l  was  computed  from  the  potential- 
independent  capacitive  "baseline"  current  between  -0.3  and  -0.7  V,  using  the  relationship 
WA  =  Ct  and  the  approximation  Cx  —  C2  (as  discussed  in  the  preceding  section).  Next, 
the  peak  potential  of  the  experimental  voltammogram,  Ep  =  -0.93  V,  referenced  to  the  Epzc 
of  Hg  in  a  non-adsorbing  electrolyte  (Epzc  =  -0.45  V  vs.  SSCE  in  NaF),  was  used  to 
compute  a  value  of  E®  for  the  lipid  monolayer  from  the  relationship  Eo  =  Ep  -  Epzc-  The 
charge  on  the  molecule  was  taken  as  +1,  and  the  parameters  T,  Ceiec>  and  z  were  given 
their  respective  experimental  values.  Finally,  the  parameters  £i/di  and  Fx  were  adjusted  to 
yield  the  best  visual  fit  of  the  theoretical  curves  to  the  data.  Values  of  £i/di  =  20  A- 1  and 
Fx  =  2.2  X  10*11  mol/cm2  were  obtained  by  this  procedure. 

The  reasonable  agreement  between  the  theoretical  and  experimental  voltammograms 
shown  in  Figure  10  supports  the  notion  that  the  peak-shaped  voltammogram  results  from  a 
conformational  transition  within  the  lipid  monolayer.  However,  as  discussed  above,  our 
assumed  structures  are  overly  simplified,  and  do  not  include  the  possibility  of  bilayer 
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formation  as  proposed  by  Nelson  and  Leermakers  [13].  The  fitted  value  of  e2/d2  =  0.241 
A-1,  corresponds  roughly  to  that  expected  based  on  a  length  of  ~16  A  [29]  and  an  effective 
dielectric  constant  of  ~3.9.  This  value  of  £2  is  in  resonable  agreement  with  that  expected 
for  the  hydrocarbon  tail  region.  On  the  other  hand,  the  shape  of  the  wave  is  not 
particularly  sensitive  to  the  chosen  value  of  ei/di,  for  values  of  ei/di  >  1,  as  evident  in 
Figure  9.  Thus,  our  fitted  value  of  Ei/di  =  20  A-1  should  not  be  interpreted  as  a  precise 
value  for  this  parameter.  However,  using  a  reasonable  value  of  di  (~2  A).  gives  a 
numerical  value  of  ei=40,  in  good  agreement  with  literature  values  [30]. 

The  value  of  Ft  =  2.2  x  10-1 1  mol/cm2  obtained  from  the  fitted  data  is 
approximately  an  order  of  magnitude  smaller  than  that  expected  based  on  molecular 
models.  Our  underestimation  of  Ft  may  indicate  that  the  charge  per  lipid  molecule  is 
considerably  less  than  the  assumed  value  of  unity  used  in  the  calculation,  or  that  only  a 
fraction  of  the  molecules  undergo  a  conformational  transition  that  results  in  the 
displacement  of  charge  within  the  interfacial  field.  As  noted  previously,  the  height  (and, 
thus,  the  area)  of  the  voltammetric  wave  has  essentially  the  same  dependence  on  Zmoi  as 
Fx*  Thus,  a  diminution  of  Zmoi  to  0.1  in  the  theoretical  curve  is  equivalent  to  increasing  Fj 
by  a  factor  of  10. 

The  inset  in  Figure  10  shows  the  interfacial  potential  distribution  corresponding  to 
the  calculated  voltammogram.  The  curves  show  that  the  transition  between  the  outwards 
and  inwards  directed  conformation  does  not  occur  until  the  potential  is  significantly 
negative  of  the  Epzc,  a  consequence  of  the  relatively  greater  stability  of  the  monolayer 
conformation  in  which  the  charged  head  groups  are  oriented  towards  the  polar  solution 
phase. 

Conclusions.  A  general  analytical  method  has  been  developed  that  allows  calculation  of 
the  voltammetric  response  of  electrochemically-inert  monolayers  that  are  susceptible  to 
field-induced  conformational  transitions.  The  model  predicts  the  presence  of  peak  shaped 
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voltammetric  features  which  are  quite  consistent  with  numerous  experimental  observations. 
The  method  provides  a  means  to  determine  the  charge  displacement  and  energetics 
associated  with  the  conformational  transition,  and  can  easily  be  adapted  to  amphiphilic 
monolayers  with  structures  and  potential-dependent  interactions  different  from  those 
employed  in  this  initial  work.  However,  a  rigorous  test  of  this  method  will  require  input  of 
potential-dependent  structures  obtained  from  direct  experimental  observations.  We  are 
currently  carrying  out  in-situ  X-ray  reflectivity  as  well  as  in-plane  surface  X-ray  diffraction 
studies  of  these  systems  which  will  allow  us  to  make  such  quantitative  comparisons  in  the 
future. 
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Table  I.  Comparison  of  Apparent  (Qapp)  and  True  (Qf)  Surface  Coverage  as 
a  Function  of  Lipid  Monolayer  Parameters. 


Qapp  (mCym2)b 


Of  (mC/m2)a 

ei/di  =  2.5  A'l 
e2/d2  =  0.14 

£l/di  =  0.5  A'l 
e2/d2  =  0.14 

ei/di  =  2.5  A-1 
e2/d2  =  0.5 

12 

11 

9 

7 

17 

16 

13 

13 

23 

21 

18 

17 

29 

27 

22 

22 

35 

32 

27 

27 

40 

38 

31 

32 

46 

43 

35 

36 

52 

49 

40 

41 

58 

54 

44 

46 

64 

59 

49 

50 

a)  Of  =  charge  corresponding  to  the  true  coverage  (zhioiFTt) 

b)  Qapp  =  charge  determined  from  the  integrated  area  under  the  calculated  voltammetric 
wave  (see  text  and  ref.  [27]). 
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Figure  Captions. 

1 .  Cyclic  voltammetric  responses  of  a  Hg/DOPC  electrode  immersed  in  a  N2-purged 
0.20  M  KCl  solution  (pH  =  8.2).  Scan  rate:  100  V/s. 

2.  Schematic  diagram  of  the  orientation  of  amphiphilic  molecules  in  response  to  the 
electrode  potential  (relative  to  the  potential  of  zero  charge,  Epzc)- 

3.  Model  system  used  to  calculate  the  charge  distribution  and  electric  field  in  the 
amphiphilic  monolayer. 

4.  Voltammetric  response  (solid  line)  of  Hg  electrode  coated  with  a  amphiphilic 
monolayer,  calculated  from  eqs.  (17)  and  (19)  using  the  following  parameters:  £« 
=  0,  Zmol  =  0.5,  Ft  =  6  X  10-11  mol/cm2,  ej/di  =  5  A'l,  e2/d2  =  0.182  A'l,  Ceiec 
=  0.10  M,  z  =  1,  £3  =  78,  and  T  =  300  K.  Fraction  of  lipid  molecules  (f  =  Fi/Tt) 
oriented  with  charged  head  group  towards  the  electrode  surface  (dashed  line). 

5.  Relationship  between  charge  density  on  metal  (Om)  and  electrode  potential  (E) 
The  dashed  line  indicates  the  shift  in  the  Epzc  due  to  the  adsorbed  monolayer 
(relative  to  a  bare  electrode).  System  parameters  used  in  calculation  are  the  same  as 
in  Fig.  4. 

6.  Potential  distribution  across  the  Hg/lipid  monolayer  interface  as  a  function  of 
applied  potential,  corresponding  to  the  voltammetric  response  shown  in  Fig.  4. 

7 .  Voltammetric  response  as  a  function  of  E°.  All  other  system  parameters  are  the 
same  as  in  Fig.  4. 

8 .  Voltammetric  response  as  a  function  of  (a)  surface  coverage,  Fj  and  (b)  electrolyte 
concentration,  Cleiec-  All  other  system  parameters  are  the  same  as  in  Fig.  4. 

9.  Voltammetric  response  as  a  function  of  (a)  ei/di  and  (b)  e2/d2.  All  other  system 
parameters  are  the  same  as  in  Fig.  4. 

10.  Comparison  of  experimental  voltammogram  (solid  circles,  data  from  Fig.  1)  and 
theoretical  response  calculated  using:  E®  =  -0.475  V,  Epzc  =  -0.450  vs  SCE,  Fj  = 
2.2  X  10-11  mol/cmA  Zmol  =  F  Qiec  =  0-10  M,  £3  =  78,  £i/di  =  20  A-1  and  (b) 
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e2/d2  =  0.241  A-1.  The  inset  shows  the  interfacial  potential  distribution 
corresponding  to  the  voltammetric  response. 
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